Production of double hypernuclei with high energy antiprotons at PANDA 
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The data available in literature, concerning the binding energy of double hypernuclei and their 
production techniques are briefly reviewed. Then, a new technique for producing double hypernuclei 
with antiprotons in flight and measuring their binding energy, proposed for the PANDA experiment 
at GSI, is investigated. Furthermore, preliminary results of the calculations for evaluating the double 
hypernuclei production and detection rates at the antiproton beam intensity foreseen at HESR are 
reported. 



1. INTRODUCTION 

The search for Double Hypernuclei (DH) (nuclei with 
two A's replacing two non-strange nucleons), started 
in the sixties with the pioneering works of Danysz et 
al. and Prowse et al, who first observed ^A^e jl| and 
A ^He |2j in emulsions exposed to K~ beams. Both ex- 
periments measured -Baa(aa Z ) = B\( A A Z) + B\( A ~^Z) 
and ABaa(aa^) = b a(aa z ) - B A { A ^Z). The measure- 
ment of these quantities is at present the only experimen- 
tal way to collect information about the A-A interaction 
and this is the reason of the stronger interest for the DH 
investigation with respect to single hypernuclei. 

Other interesting features of DH are: a) their nu- 
clear structure, with two levels both filled by strange 
hadrons; b) the possibility to explore the existence of the 
H dibaryon (a double strange system of 6 quarks pre- 
dicted several years ago by Jaffe 3] ) ; c) the information 
about the levels of the S~-atom, that is formed during 
the DH production process. In spite of the interest of 
these topics, only few experiments have been devoted 
to DH so far, mainly because of the difficulty to pro- 
duce them. The traditional way is to produce a double 
strange S _ particle via the strangeness exchange reac- 
tion K~(p, 'Er)K + between a K~ meson and a proton 
bound in a nucleus, giving K + and S~; this last, if de- 
celerated to rest before decaying, can be captured inside 
a nucleus in which it interacts with a proton eventually 
releasing two A's in two hypernuclear levels. The 2 _ en- 
ergy may be of order of some hundred MeVs, thus the 
slowing down process takes long time in ordinary matter. 
This explains why the probability of DH formation is so 
low: in fact, Danysz and Prowse Q observed one DH 
event each, among a K~ number of around a million. 

In more recent experiments at BNL-AGS (E885) and 
KEK (PS E176) about 2 • 10 4 and 800 stopped S"'s re- 
spectively were observed, as reviewed by Pochodzalla Q . 
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Also at the planned Japanese Hadron Facility (JHF) the 
number of produced S _ 's is expected to be of order of 
some thousands. 

Furthermore the results obtained from experiments in 
terms of the above mentioned binding energies are too 
much different, clearly outside the experimental errors 
(see again |4j). Thus, while the existence of DH seems 
ascertained, an increase of statistics in the S~ produc- 
tion is mandatory in order to step forward in their un- 
derstanding. 

Even though the traditional double strangeness ex- 
change reaction will be pursued in the next future at JHF, 
as already mentioned, nevertheless the advent of new an- 
tiproton facilities, like HESR and JHF (second phase), al- 
lows the scientific community to explore alternative tech- 
niques of S~ production from antiproton beams. A first 
one was proposed by Kilian [5| who suggested the annihi- 
lation of antiprotons at rest to produce K~ : s which sub- 
sequently could exchange strangeness with the residual 
nucleus, giving a S - hyperon. This technique requires 
a low energy antiproton machine and could be imple- 
mented, for instance, at CERN-AD. 

But also in high energy antiproton machines like HESR 
and JHF, the S~ hyperon may be directly produced in 
the antiproton annihilation reaction together with its an- 
tiparticle. The advantages are: 

1. an higher rate with respect to the double reaction in- 
volving the annihilation at rest; 

2. the presence of the antiparticle, to be used for tagging 
purposes. 

The only disadvantage is the high momentum of the 
S~'s which need to be strongly decelerated before being 
captured at rest. 

In this work we present the first preliminary results 
of a Monte Carlo simulation based on a simplified Intra- 
Nuclear Cascade Model, and performed to explore the 
rates of the produced and stopped S - 's and their proba- 
bility to be captured forming a DH, as may be expected 
in the future PANDA experiment at the HESR-GSI fa- 
cility. 
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In sec.|2 the 3~ production techniques are illustrated; 
in sec. OH the details of the antiproton annihilation into 
the 3~'s are discussed; in sec. 0] the slowing down pro- 
cesses of S~ inside the residual nucleus, in the target and 
in a second target devoted to the DH formation are anal- 
ysed; in sec. the results of our simulation in terms of 
DH rates are reported, and finally in sec. the conclu- 
sions are drawn. 

2. DOUBLE HYPERNUCLEI PRODUCTION 
WITH ANTIPROTONS 

The simplest way for producing double strange nuclear 
systems is the DH direct formation through K~K + re- 
actions; on the contrary, an indirect method is the for- 
mation using 3 _ capture at rest. The latter is known to 
have a higher rate jg]. 

A classical technique, adopted in the El 76 experiment 
at KEK [7J,|8|, is based on a K~ kaon beam that induces 
the quasi-free nuclear reaction K~ + p — > K + + S~, thus 
producing a relatively high number of fast 3~ hyperons. 
These hyperons subsequently slow down in an emulsion 
stack and some S = — 2 hypernuclei are formed following 
the S _ capture at rest. In such kind of experiments, the 
emulsion is also the detector, with high enough resolution 
for tracking the formed hypernuclei and the products of 
their decay as well. A further improvement has been ob- 
tained in E373 at KEK |6j, through a K~ beam (with mo- 
mentum ~ 1.66GeV/c) impinging on a diamond target 
in which the same reaction previously sketched occurs, 
though at far higher rate: in the experimental setup, the 
diamond block was followed by an emulsion stack as for 
E176 and the tracking detector was a fiber-bundle sys- 
tem. In all these experiments, the S~ hyperon is cap- 
tured in a light nucleus (i.e. carbon, nitrogen or oxygen) 
of the emulsion and therefore the observed double hyper- 
nuclei are A^He, AA.Be and /$J5. 

Let us consider now the new technique proposed for the 
planned PANDA experiment at HESR-GSI; it relies on 
the same method of stopped S _ 's, the main difference 
being that the hyperon production reaction is realized 
through an intense beam of antiprotons on a target of 
A X^ nuclei (where A is the mass number, Z is the atomic 
number and N is the neutron number) . The two reactions 
of interest are 

p + A x£^ 3° + H- + A - 1 X^ 1 + mesons , (1) 

and 

J 3 + A x£^S+ + S-+ A - 1 X'z_ 1 + mesons. (2) 

In reaction^ the energetic antiproton annihilates on a 
neutron bound inside the X nucleus, while in reaction [3 
it annihilates on a proton. 

The fast 3 _ slightly slows down in the production tar- 
get itself, but the main part of its energy-loss process 



occurs in a second separate target, that is also the place 
in which a statistical collection of double hypernuclei is 
supposed to be formed with some probability. In both 
reactions n and El beside the S = — 2 hyperon S~, an 
anti- hyperon (3° or S + ) is released too, according to the 
strangeness and baryon conservation laws of strong in- 
teractions: this antiparticle can play a crucial role in the 
experimental detection of the whole double hypernuclei 
formation, as it may be used for trigger purposes. 

Before proceeding, we want to point out that a still 
different technique has also been proposed Q, in which 
an antiproton interacts at rest with the A X^ nucleus. A 
possible reaction is 

p + A X£-> K*- +K° + A - 1 X^- 1 + mesons, (3) 

where K*~ is the K~ resonance at 892 MeV, that subse- 
quently interacts with the residual nucleus, giving 

K*- + A - 1 X^- 1 -> 3- + K° + A - 2 X£- 2 + mesons , (4) 

where the 5~ hyperon emerges with a momentum be- 
tween 250 and 800MeV/c approximately. In comparison 
with the method based on reactions ^ an d H this tech- 
nique has the disadvantage that the S~ production oc- 
curs in two distinct steps (reactions [3] and , each one 
characterized by a somewhat low probability. 

3. H PRODUCTION 

Let us start considering the S~ production from re- 
action ^ A highly energetic antiproton p interacts with 
a neutron n bound in a nuclear potential well. In what 
follows, we assume that the neutron be free; indeed, we 
expect that the total antiproton energy be of order of a 
few GeVs, while the average binding energy of a nucleon 
is roughly 8.8 MeV for the most tightly bound nuclei (i.e. 
56 Fe, 58 Fe and 62 Ni). 

If we consider p momenta below tt production thresh- 
old, the previous discussion allows us to consider the re- 
action as a two-body process, with the manifest advan- 
tage that its simple kinematics can be treated analyt- 
ically, starting from the relativistic four-vector energy- 
momentum conservation; a similar discussion also applies 
to reaction Hence, the two reactions simplify as 

p + n^3°+S-, (5) 

and 

p + p^E++E-, (6) 

and from now on, we shall refer to them both as 
Strangeness Creation Reactions, namely SCRs. 

The calculated threshold for the two SCRs (neglecting 
the slight difference between the proton mass and the 
neutron mass) is located at an antiproton momentum of 
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Figure 1: Plot of the hyperon momentum P(S~) against 
the exit angle 9(S~) after SCR, considering three different 
values of the antiproton momentum, namely Po = P thr = 
2.62GeV/c, P = 2.81GeV/c and P = P* hr = 3GeV/c. 
Both first and second kinematical solution are shown together. 
The curve labelled Po = 3 GeV/c has been adopted in our cal- 
culations. 



pthr _ 2 02 GeV/c. Furthermore, the upper limit for the 
two-body exit channel (i.e. below pion production) is 
P^ hr = 3GeV/c. 

Let be the angle between the direction of the 

incident antiproton p and the direction of the outgoing 
S~ hyperon, in the laboratory frame of reference. As 
it is well known from the relativistic kinematics of two- 
body reactions with threshold, the 0(E~) angle ranges 
between and a maximum, 9 max , that is a function of the 
antiproton momentum P . Furthermore, for any given 
value of t9(S _ ), two different momenta of the S~ hyperon 
are allowed, whose relative probability strongly depends 
on the dynamics of the reaction itself. All these peculiar 
characteristics are shown in fig. referring to the SCR 
of eq. El the dependence of the maximum angle # max on 
Po looks evident. At the antiproton momentum Po = 
3GeV/c, we obtain that the maximum angle is # ma x = 
0.3 rad. 

The argument that led us to the choice of Po relies 
on the results of the quark-gluon string model [t| which 
predicts that the p + p — * S+ + S~ reaction cross sec- 
tion a shows a maximum at an antiproton momentum 
of around 3 GeV/c; in this case, the estimated maximum 
value is a sen. ~ 2 /ib. Then, in our following calcula- 
tions, we adopt Po = 3 GeV/c and we suppose that the 
same total nuclear cross section, ctscr = 2 /ib, does hold 



for both SCRs of eqs. |S] and HO Lying on these physi- 
cal hypotheses we obtain that, as far as the first (most 
probable) kinematical solution is concerned, the S~ mo- 
mentum ranges from 2129MeV/c in the forward direc- 
tion, to 1301MeV/c at the limit angle while, considering 
the second solution, it spans from 877MeV/c to about 
1301 MeV/c respectively, as it can be seen in fig. ^ 

Because of the high average momentum, the energy- 
loss process is of central importance in the planned ex- 
periment, and therefore the evaluation of the S~ slowing 
down is worth while; we have performed it through a 
Monte Carlo simulation, as described in sec. As a 
starting point, we generate 2 • 10 5 S - particles produced 
from SCR for each solution, assuming that reaction[S]oc- 
curs through a S-wave in the centre of mass frame of refer- 
ence, namely assuming an uniform spherical distribution. 
According to the 3 GeV/c curve in fig-Q the momentum 
spectra of the simulated S - hyperons are presented in 
fig. |3 for both kinematical solutions. 

For tagging purposes, the momentum spectrum of the 
anti-hyperon is of great interest; our computer simulated 
results are reported in fig- El The S° (or S + ) momentum 
ranges from 871MeV/c to 1798MeV/c, with an average 
emission angle of 0.262 rad, in first solution, and from 
1799MeV/c to 2122MeV/c in second solution, with an 
average angle equal to 0.125 rad. Furthermore, the max- 
imum angle is 0.3 rad; therefore, the tagged anti-hyperon 
is expected to be released very close to the antiproton 
beam direction. 



4. 3 - SLOWING DOWN 

Let us start by noticing that the S~ slowing down 
process consists of two distinct steps: 

1. a sequence of nuclear elastic scattering events with 
some of the A — 1 nucleons of the residual nucleus in 
which the annihilation has occurred; 

2. the energy loss by ionization in the ordinary matter of 
the two targets. 

We shall discuss each step separately. 

4.1. Hyperon-nucleon elastic scattering 

We assume, on the basis of ref. 0], that the total 
elastic cross section for the and 'ETn scattering pro- 
cesses be <je ~ 10 nib. Moreover, we assume a differ- 
ential elastic cross section d<7E/dfi oc exp(P ■ t), where 
t is the second Mandelstam variable, and B is suitably 
taken as 5GeV -2 0]. Furthermore, we model the nu- 
cleus as a homogeneous sphere of nucleons of radius 
P nuc = R ■ (A - 1)V3 ( Ro a 1.35 fm). 
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Figure 2: Simulated P(H ) momentum spectrum after SCR, for both first and second kinematical solution. We assume an 
isotropic angular distribution in the centre of mass frame of reference. 
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Figure 3: Simulated P(H°) momentum spectrum after SCR, for both first and second kinematical solution. 



We perform the numerical simulation of this first slow- 
ing down step in the framework of an INC-like (Intra 
Nuclear Cascade) model jj^i starting from the distribu- 
tion of fast hyperons shown in fig. The basic hypoth- 
esis of this model is that the (A — 1) residual nucleus 
does survive for a time longer than the time spent by 



the hyperon during its intra-nuclear path. Furthermore, 
the scattering exit angle of S~ after each nuclear scat- 
tering is chosen uniformly in the centre of mass frame of 
reference. 

As far as nuclear scattering events are concerned, we 
treat these processes as instantaneous since the calcu- 
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lated average spent time is of order 10 22 s (proper time), 
while the S~ mean life, in its rest frame of reference, is 
t ml = 1.639 -10- 10 s. 

We obtain that, from the 2 • 10 5 simulated random 
walks, 9.1% of the hyperons scatter in the residual nu- 
cleus at least once if we consider the first kinematical 
solution (the maximum number of consecutive scattering 
events being 3), and 31% considering the second solu- 
tion (and the maximum number of scattering events is, 
correspondingly, 8). 

The main physical effects due to hyperon-nucleon 
elastic scattering are twofold: on one hand, the 
P(!H _ ) momentum spectrum is modified with respect 
to that of fig. showing a low momentum tail from 
1301MeV/c down to ~ 96.6MeV/c (first solution) and 
from 877MeV/c to ~ 18.7MeV/c (second solution); on 
the other hand, the angle 9 nuc between the exit direc- 
tion and the antiproton direction (8 nuc = if the 
hyperon does not scatter at all), after the nuclear scat- 
tering process, spans from up to 1.48 rad and 2.38 rad, 
respectively, which are much higher than the limit an- 
gle # m ax defined in sec. [31 Thus, this first step provides 
a non negligible fraction of hyperons with relatively low 
momentum and large angles. The next step, involving 
the passage of particles through matter and the ioniza- 
tion energy-loss, further modifies the momentum spec- 
trum, without modifying the 6* nuc angle. Due to the low 
probability, the nuclear scattering of the hyperon inside 
the target has been neglected. Thus, the final flight in 
the experimental apparatus is a simple straight line. 



4.2. Energy-loss by ionization 

In the study of the energy-loss during the passage of 
particles through matter, the geometry of the experimen- 
tal setup plays a crucial role, because it determines the 
effective hyperon path as a function of the exit angle 9 nuc , 
and thus the slowing down effectiveness. 

We adopt here a simple geometry consisting of two 
distinct targets: the first one, in which the SCR takes 
place, is a thin parallelepipedal wire of Gallium or Gold 
(4 cm long and 5 /im x 5 /jm square) orthogonal to the 
beam direction, while the second one is a diamond paral- 
lelepiped (4 cm x 4 cm as transverse dimensions and 2 cm 
thick) with a hole of radius 0.25 cm in its centre. The 
two targets are separated by a 1 mm vacuum gap. 

The physics of our simulation essentially relies on the 
Bethe-Bloch equation, which expresses the mean rate of 
energy-loss per unit length: we take it from |l3| and we 
report it here for the sake of discussion 
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where Ek is the hyperon kinetic energy, p is the density 
of target material, K = 0.3 MeV • cm 2 , Z and A are the 
atomic number and the mass number of medium, [3 is 
the dimensionless velocity of hyperons, m c is the electron 
mass, 7 = (1 — /3 2 ) -1 / 2 is the relativistic factor, / is the 
mean excitation energy, Mq is the hyperon rest mass and 
c is the speed of light. The parameters C and 6, usually 
negligible at high energy, play an important role in our 
similation because the final stage of the slowing down 
process takes a major fraction of the total time spent by 
the hyperon before being stopped; the value of 5 has been 
taken from the accurate work of Sternheimer et al. [lif 
and reads 



2Xlnl0 + C 
.) = ( 2X\nlO + a(X 1 - X) 7 




X > Xi 
c X < X < x 1 
x<x 



(9) 

with X = ln/?7, and where C, a, m, Xq, Xi are taken 
from ref. for the diamond target. 

The total time T elapsed during the energy-loss pro- 
cess is comparable with the S~ mean life, tml- thus, a 
fraction of hyperons decays prior to complete stopping. 

In order to evaluate the number of surviving S — 's, the 
following procedure has been adopted. Each S~ has been 
assigned a lifetime r according to the exp(— t/tml) dis- 
tribution, in the frame in which the particle is at rest; at 
each interval of lost kinetic energy dEx , the correspond- 
ing time interval dt (in the laboratory frame) spent by 
S~ to travel a path dx has been transformed into the rest 
frame in order to obtain the total proper time elapsed T, 
and compare it with the proper lifetime r. Thus, eq. [7| 
should be solved provided that the condition 



T= / dty/l ~f3 2 <T, 



(10) 



be satisfied, where is the final time after complete 
stopping of the hyperon in the laboratory frame of refer- 
ence. 

Eq. can be formally expressed through 
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Furthermore, the hyperon momentum P depends over 
Ek as 



P{E K ) = y/E K {E K + 2M„) . (12) 

From eas. 1101 ITTl and IT2l we straightforwardly obtain 

,o 1 

R=- dE K —— (13) 
J E i s g{E K ) 
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JUT pO i 

T = / dE K - — < t , (14) 

c J E u K g(E K )P{E K ) ~ K > 

where 7? is the stopping range and is the hyperon's 
initial kinetic energy. Performing a numerical integration 
of eq. 1131 at each step the corresponding time integral 
of eq. El is calculated and checked whether T is lower 
or greater than the lifetime r. In the former case the 
integration proceeds until complete stopping, otherwise 
the hyperon decays along its path. 

5. SIMULATION RESULTS 

The simulation procedure described in the previous 
section has been applied to the above mentioned geom- 
etry (a wire as S~ production target followed by a sec- 
ond one as hypernuclear target); furthermore we run the 
simulation twice, considering two different materials (gal- 
lium and gold) for the wire, in order to check the effec- 
tiveness of the atomic weight on hyperon production. For 
both materials, 2 • 10 5 S~ 's have been generated in a spot 
of radius r = 10 /im on the surface of the wire. 

The calculated parameters of our simulation are: 

• the fraction / s top of hyperons stopped in the second 
target; 

• the fraction /d ec of hyperons decayed before stopping; 

• the fraction /i oss of hyperons lost in the gap between 
the two targets or in the central hole of the second 
target. 

Since reactionnjcan produce S~'s of two different mo- 
menta, both kinematical solutions have been simulated 
separately. 

The most relevant parameters for the experiment are 
the total number N stop of stopped and detected S~'s 
per second and the number AAa of produced and de- 
tected double hypernuclei per second. In order to cal- 
culate these values from the simulated fractions already 
defined, one has to make a few reasonable assumptions 
about the luminosity of the machine and the detection ef- 
ficiency; moreover, other physical parameters mentioned 
below have been roughly estimated just to get a first in- 
sight of the rates, though further accurate evaluation is 
needed. 

As already discussed in sec. 03 we have assumed a total 
cross section <tscr ~ 2 fib for the SCR reaction between 
an antiproton and a free nucleon at 3 GeV/c; accordingly, 
the corresponding reaction cross section of reaction ^ oc- 
curring on a neutron inside a nucleus reads 

Zscr « vscr ■ A 2 ' 3 ■ , (15) 



Table I: Results of the Monte Carlo simulation of the double 
hypernuclei production, considering two materials (gallium 
and gold) for the hyperon production wire. Each kinematical 
solution has been investigated separately, and the correspond- 
ing results are shown. (See text for the physical meaning of 
the different parameters). 
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where the surface annihilation scaling has been taken into 
account [lij . 

Furthermore, we assume a luminosity of the antipro- 
ton beam of roughly C w 10 32 cm _2 s _1 , as realistically 
expected at HESR, and a reconstruction efficiency for the 
whole detector ek ~ 0.5. 

Following refs. 0, ^| , we also assume that the con- 
version probability (i.e. the probability for the hyper- 
onic S~-nucleus conversion into a AA-hypernucleus) per 
stopped S~ be paa ~ 0.05, and that the level population 
fraction be plp ~ 0.1; the probability of transition per 
event has been taken as pt ~ 0.5 and the 7-photo peak 
efficiency as e 1 w 0.1. 

Thus, the number of stopped S~'s and of formed AA- 
hypernuclei can be calculated through 

Astop = £ ■ ^SCR ■ /stop • £k , (16) 

and 

AAa = A^top • paa • Pt ■ Plp ■ e 7 ■ (17) 

The final results are reported in tab.[I]fOT both gallium 
and gold wires and both kinematical solutions. 

The relative weight of each solution is at present totally 
unknown as no measurements of the differential cross 
section for reactions and El exist nor theoretical esti- 
mates as well. Hence, the actual value of AAa is located 
between the two values of first and second solution re- 
ported in tab. [I] Being conservative, one may take the 
value of the first solution only; nevertheless, even in this 
pessimistic hypothesis, the total number of stopped S~'s 
is around 1.19 • 10 6 per month and the total number of 
formed double hypernuclei is around 300 per month con- 
sidering the gallium wire, while we obtain 3 • 10 6 and 750 
respectively when considering the gold wire. 

On the other hand, if the p annihilation occurs through 
reactionElinstead of reaction^] the processes of H~ slow- 
ing down and double hypernuclei formation are in prac- 
tice identical to those considered so far: therefore, it is 
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reasonable to expect that the total 5 production and 
double hypernuclei formation should be enhanced by a 
factor A/ (A - Z) ~ 1.8. 

As a first indication from the values of the stopped 
E~'s, one can observe that the weight of the nucleus 
strongly influence the hypcron slowing down through nu- 
clear scatterings. 

Of course these results depends on the estimates of 
the machine, detector and physical parameters chosen for 
our simulation, but the values assumed here seem quite 
realistic indeed. In fact, HESR design's first priority is 
the high luminosity and the PANDA apparatus will be 
an ensemble of several detectors with good reconstruction 
capabilities. 

6. CONCLUSIONS 

A preliminary evaluation of the rates for the complete 
stopping of a S~ hyperon, produced through a high en- 
ergy antiproton interacting with a gallium or a gold nu- 
cleus, in a diamond (carbon) target has been performed 
using a Monte Carlo technique based on an INC-like 
model. The chosen geometry is suitable for the PANDA 
experiment project at the future HESR machine at GSI. 

The results, relying on the expected performances of 
the apparatus, show the feasibility of a large produc- 
tion of double hypernuclei: our preliminary estimates, 
in terms of stopped S _ 's and formed double hypernuclei, 
are larger than the previous data existing in literature 
and also larger than the expectation of other future ex- 
periments and machines. 

A lot of calculations are of course still necessary in 
order to fully design the geometrical arrangement, the 
sizes and the materials of the targets, in order to optimize 
the production and to meet the beam requirement of the 



machine, and this will be the future work for the next 
months. 
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